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The current-induced motion of a domain wall in a semicircle nanowire with applied Zeeman field 
is investigated. Starting from a micromagnetic model we derive an analytical solution which char- 
acterizes the domain-wall motion as a harmonic oscillation. This solution relates the micromagnetic 
material parameters with the dynamical characteristics of a harmonic oscillator, i.e., domain- wall 
mass, resonance frequency, damping constant, and force acting on the wall. For wires with strong 
curvature the dipole moment of the wall as well as its geometry influence the eigenmodes of the 
oscillator. Based on these results we suggest experiments for the determination of material param- 
eters which otherwise are difficult to access. Numerical calculations confirm our analytical solution 
and show its limitations. 

PACS numbers: 75.60.Ch, 72.25.Ba, 76.50.+g 



I. INTRODUCTION 

Field-driven dynamics of magnetic domain walls have 
been intensely studied over the last decadespi*^ The 
topic has recently regained interest by the discovery 
that spin-polarized currents of high density can al- 
ter magnetization configurationaSiii^ and move domain 
.^g^jjg|7j8j9ji0jii Current-induced magnetic switching is 
viewed as a promising solution for the realization of 
magnetic random access mcmoriesSiiSii^, while current- 
induced domain-wall motion has potential application 
in spintronic data storage devices, e.g. in the race- 
track memorjii^ or data transfer schemes li^iiSiii Sev- 
eral models of current-driven magnetization dynamics 
have been established to explain the electronic origin 
of current-induced magnetization changes and to predict 
their effects i^i^iiSiiSiSfliSi At first it was assumed that for 
finite domain walls the spins of the conduction electrons 
adiabatically follow the local magnetic momentsiiSiSS 
Later the theoretical model was extended to include a 
non-adiabatic mismatch between the current polarization 
and the direction of magnetizationiiSiSfliSi 

The measured and calculated velocities of current- 
driven magnetic domain walls in thin nanowires vary 
by several orders of magnitude even for the same 
material iSiiSSiSliS^ While it has been originally sug- 
gested that the discrepancy could be due to thermal 
activationii42SiSS42L2& or surface roughness^^, it has re- 
cently been found that the domain-wall velocity depends 
on the type of the domain wallSSiSS which can be changed 
by a spin-polarized current ASSiSSiSi Recently it has been 
observed that the velocity of field-driven domain-wall 
motionS^ can be altered by ± 100 m/s by a pulsed spin- 
polarized current^ and that the motion can even be 
halted completelyi^ It is now assumed that the adia- 
batic term is largely responsible for the acceleration of 
the domain wall while the non-adiabatic term will cause 



the wall to continually move?^ It has been shown that 
domain-wall oscillations excited with an ac current at its 
resonance frequency require current densities one to two 
orders of magnitude less (10^" A/m^, see RefA^^) than 
for pulsed excitations (lO^-lO^^ A/m^, see RefALii*2&). 

Here we show that the harmonic-oscillator model fol- 
lows naturally from a micromagnetic model that de- 
scribes the excitation of transverse walls in thin narrow 
j.jjjgg|36j37j38 Solving analytically the Landau-Lifshitz- 
Gilbert equation extended by the current corrections due 
to Zhang and LiSi we are able to express the proper- 
ties of the driven oscillator by the quantities determining 
the micromagnetic model. A comparison of the numeri- 
cal calculations with our analytical solution confirms the 
importance of the geometry due to the curved wires. Fi- 
nally, we suggest experiments which can determine the 
values of the non-adiabatic spin torque and the Gilbert 
damping parameter. 



II. MODEL 

Figure ^ shows a ferromagnetic semicircle nanowire 
with a domain wall at its bottom placed in an exter- 
nal magnetic fieldi^ The wall is excited by an oscillating 
current flowing between the two contactsiS 

The magnetization dynamics of a magnetic wire is 
well described by the Landau-Lifshitz-Gilbert (LLG) 
equationiiS In the presence of a spin-polarized current 
density j, the interaction between the itinerant elec- 
trons and the magnetization M leads to an extension 
of the LLG equation. This extension was derived from 
a quantum mechanical model by Zhang and Li^ Their 
semiclassical approximation results in the extended LLG 
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FIG. 1: Scheme of the semicircle nanowire with radius r in 
a magnetic field H. The static magnetization in the absence 
of a current is indicated by small arrows. The two rectangles 
under angles /? = ±45° are the electrical contacts. 



equation (in Gilbert form) 
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with the gyromagnctic ratio 7. the Gilbert damping pa- 
rameter a, the saturation magnetization Af^, and the ra- 
tio between exchange relaxation time and spin-flip relax- 
ation time ^ = Tex/^gf. The effective magnetic field H^q 
includes the external as well as the internal fields. In this 
model the spin current is sensitive to the spatial inhomo- 
geneities^of the magnetization with a coupling constant 
where P denotes the spin polarization of 



"0 ~ eM,(l+i^, 

the current and ^lb is the Bohr magneton. 

Since the saturation magnetization is constant for a 
given material at fixed temperature, Af is perpendicular 
to ^ and Eq. d 

can be reformulated to an explicit 
equation of motion for the magnetization 
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This 



with the abbreviations 7' = 

equation is the starting point for the analytical as well as 
for the numerical calculations presented in the following. 



III. ANALYTICAL CALCULATIONS OF THE 
STRAIGHT WIRE 

For the analytical treatment of Eq. ||2Jl we transform 
the semicircle wire in a homogeneous Zceman field to 



a straight wire in a spatially variing field. The wire is 
directed along the a;-axis. The direction of the mag- 
netization will be expressed in a polar spin basis M ~ 
Afs(cos 0, sin0 cos 0, sin0 sin (^). In the absence of elec- 
tric current and external magnetic field the energy of a 
domain wall within the wire is 



E = S 



mix) 
dx 



K sin^ I 



dx. 



(3) 



where denotes the angle between the wire axis and the 
magnetization. A and K denote the exchange and the 
shape anisotropy constant. This functional can be mini- 
mized by the well known Neel wall described by the angle 



= TT — 2 arctan 



(4) 



The center of the wall is at positi on X and the width 
of the domain wall is A = ^ A/K. From Eq. (0J two 
expressions 



cos 9 = tanh 



x-X 
X 



1 



smf = 



cosh I 



(5) 



can be derived which will be useful in our further calcu- 
lations. 

In the presence of an external field iJ^xt ^^'^ demagne- 
tization energy K± sin^ 9 sin^ (j) caused by the rotation of 
the wall around the wire axis can no longer be neglected. 
We include the external field perpendicular to the wire 
into the shape anisotropy K^. The energy functional in 
Eq. Q has to be extended to 



E = 



KsmH + Ai^ 
ox 



dV 



j \Ki_ sin^ 9 sin^ (j) — /xoA'fs-ffext(^) ^] 



(6) 



Here we have restricted ourselves to an external field par- 
allel to the wire. Also the crystalline anisotropy has been 
neglected^i From the energy functional in Eq. © we 
derive the effective magnetic field through the relation 
H = — — — 

tiff Mo SM ■ 

From the extended LLG equation ^ in the polar spin 
basis we obtain 



7' 



5E 
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+ sin(6i)6' (1 + OJ • V0 - 6' - a)] ■ V9. 



(8) 



Assuming that the moving wall stays a Nccl wall we 
can describe its motion, following the description of 
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Schreyer and Walker—, by two dynamical variables: the 
position of its center X and its angle around the wire axis 
(l){x) = (f) that is uniform along the wire. With Eq. © 
and Eq. © we get from Eq. Q and Eq. © 



sin(6l) 



X 



ai sin(6')iJext(a; 



sin(6') sin(0) cos(0) 



6' 
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and 
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— 2 sin(0)7'ai^x sin((/)) cos((/)) 
+ ^ sin(6l) cos(6i) sin2(0). 
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(10) 



Note, that X and </> depend on the position a; along the 
wire. In the following wc show that a solution consistent 
with our initial assumptions exists for small exitationsj^ 
Note, that this condition holds for realistic current den- 
sities. 

Assuming that -ffext(^) varies slowly on the length 
scale of the domain- wall width A = A/ K, sin 9 is re- 
placed by a (5- function 'k\5{x — X) in view of Eq. ||SJ). 
Also we neglect terms which are nonlinear in cj). This ap- 
proximation holds for angles </) smaller than about 10°. 

The equations of motion for the domain wall then be- 
come 



X = -X2j'Ka 



1 



^qMs 



A7'ai/ext(^) 



fe;-(l 



and 



= 7'^fcxtW-27'a/^±0 



A 



(11) 

(12) 



These equations are general equations of motion with a 
time dependent current density j. 

In the limit of a steady current and a homogeneous mag- 
netic field one can calculate the initial velocity of the wall 
by setting = 0, the initial condition of the Ncel wall. 
This leads to the initial velocity Xi ~ — A7'aiJgxt~^j (1+ 
a^)j which is exactly the value obtained by Zhang and 
LiiSi The terminal velocity Xf = — (A7i7Qxt + ^jO) /o; 
is calculated by setting = 0, i.e., stationary motion. 
This velocity is also identical to the one calculated by 
Zhang and Li. Similar relations have recently been found 
by Dugaev et ali^ 

The domain- wall mass is obtained by determining sin 
in the absence of electric currents and external fields. 
From Eq. (|11|) we obtain for the stationary motion 



X 



(13) 



Inserting this result into the 0-dependent part of the 
domain-wall energy in Eq. © and comparing with the 
energy E of the domain- wall quasiparticle 



-mA:' 
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15^ 

2 xy^Kj_ 



X' 



we arrive at the domain- wall mass 



(15) 



Note that this result relates the phenomenological 
domain-wall mass of a tail-to-tail Neel wall to the mi- 
cromagnetic material parameters. 

In the case of a curved wire the projection of a uni- 
form external field along the wire is given by i?ext(^) ~ 
-ffosin(x/r). Transferring this to our straight wire 
model, at small displacements of the domain wall (A << 
r) the wall is exposed to the external field i?ext ~ 
HoX/r. Then the equations of motion become a sys- 
tem of two coupled linear differential equations of first 
order: 



X 



-XaHa- 



-X2K^-^ 
-2a A" 1 



A 



-{1 + aO 

A 



(16) 



Except for the non-vanishing first matrix element 
—XaHo/r these equations are equivalent to those of a 
driven harmonic oscillator. For a time dependent cur- 
rent density of the form joe*^* the general solution 



fX{t) 



A, 



A_ 
0_ 



-Ft—iujft 



1 



F 



(17) 



consists of an exponentially damped starting configura- 
tion with the initial conditions described by X± and (j)± 
and a current-driven oscillation with the driving force F. 
The damping constant 



XHo 
2r 



(18) 



depends on the ratio of applied magnetic field and ring 
radius. It represents the restoring force acting on the 
domain wall. This dependence of the damping constant 
on the restoring force expresses that the damping is spa- 
tially dependent. This also leads to a second term in the 
frequency of the free oscillation 



l2Y^HoXKj_ 



Affo 
2r 



(19) 
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FIG. 2: (Color online) Schematic illustration of the magneti- 
zation in the Neel wall (solid red arrows) in a straight wire of 
width w and thickness t. H^^ and Ha are the parallel compo- 
nent of the external field and the anisotropy field, respectively. 



which is different from T. Hence the resonance frequency 



UJr 



r2 = 



-(l + a2) 



(20) 



depends explicitly on the Gilbert damping a and differs 
from the resonance frequency of a normal harmonic os- 
cillator 




uja = \ — = 



fioMsr 



(21) 



by the factor \/l + o? . The constant D is given by £> = 
Fh/X where Fh is the force on the domain wall due to 
the external magnetic field. The force 



F 



mbjjoe 



* r 1 



1 + a- 



iH, ] ex 



(22) 



torque the phase at resonance frequency between the cur- 
rent and the magnetization in z-direction is 90° when the 
ratio ^ of exchange and spin-flip relaxation time equals 
the Gilbert damping parameter (^ = a). The phase can 
be used to find out whether a non-adiabatic spin torque 
exists and to determine the value of ^ in comparison to 
the damping parameter a. 

The influence of the adiabatic torque on the position 
of the wall is obtained by setting ^ = in Eq. (|22|l . The 
x-componcnt of the force F due to the adiabatic torque 
is proportional to the time derivative of the current den- 
sity. Therefore, the adiabatic torque does not accelerate 
the wall when the current does not change in time. This 
explains the observation of Zhang and Li2i that with- 
out a non-adiabatic spin torque a domain wall subjected 
to a steady current stops moving. In contrast the non- 
adiabatic contributions to the force are proportional to 
the current density as well as to its derivative. 

In Eq. p7(l the starting configuration depends on <j)± 
and X±. The equation that follows from decoupling of 
Eq. (Unjl 



a 
2A 



4rK, 



±i- 



2XYK_i_ 



(23) 



connects (j)± with X±. Hence we have two parameters 
left for our starting configuration as expected for an os- 
cillation. 

With the above analytical model we are able to derive 
the hitherto phenomenological oscillator model^ and to 
express its characteristics by the micromagnetic material 
parameters. Likewise, the measurement of the domain- 
wall motion allows the determination of micromagnetic 
quantities. 



IV. CURVED WIRES 



induced by the current depends on the frequency Q of 
the applied current. 

The terms in Eq. H22|l can be understood as direct 
forces due to the spin torque and the precessions of 
the magnetization in the external and anisotropy field 
depicted in Fig. |21 The terms proportional to ifi ex- 
press the current-induced spin torque. They are the time 
derivatives of the inhomogenities in Eq. (|16|l . The Hq- 
dependent term is a result of the precession of the mag- 
netization in the external field which causes a rotation 
of the wall around the wire axis. The precession in the 
anisotropy field, described by the K±-term in Eq. (|22|1 . 
causes a change of the wall velocity. 

Except that the force depends on the frequency fl of 
the applied current the result for the domain wall dis- 
placement is equal to the one in a harmonic oscillator. 
With increasing the force increases and its phase shifts 
up to 90°. In the absence of a non-adiabatic spin torque 
(^ = 0) current and domain-wall displacement at res- 
onance have opposite sign. In case of a non-adiabatic 



For curved wires in a homogeneous magnetic field its 
component perpendicular to the wire has to be taken into 
account. Also the change of the magnetization due to the 
curvature becomes important. To include the perpendic- 
ular field we calculate the force on the domain wall as 
the spatial derivative of its Zeeman energy. The total 
magnetic moments parallel to the wire 

m(| = J Ms cos[e{x)]dV = -2MsSX (24) 

and perpendicular to the wire 

m^ = j MsSm\6{x)]dV = ttMsSX (25) 

are volume integrals over its magnetization that are read- 
ily calculated using the relations in Eq. (O. Note that 
TO|| is the magnetic moment of an abrupt domain wall. 
With the magnetic field Hq in y-direction the Zeeman 
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energy can be written as 



given by 



sin(/?)d/3 - j sin(/3)d/3 
L-2 P« 



dr' 



(26) 

where /3o = 7- is the angle of the position of the domain 
wah (see Fig. ^ and r, w, and t are the radius, the width, 
and the thickness of the wire. We get 

E\\ = -2fioMsSrH cos{(3o) = 2fioM,SHY (27) 

with the cross section S = wt. One recognizes that the 
energy is equivalent to the energy of a monopole with 
magnetic charge Qm = 2^qMsS. For small domain- wall 
displacements we can write the cosine in Eq. (|27|l as a 
Taylor series up to second order in X and get 



-2fioMsSrH 1 



2r2 



(28) 



The monopole has been included in the above calcula- 
tions as well as in the calculations of Saitoh et aliS The 
perpendicular magnetization contributes to the Zecman 
energy 



Ej 



-fiQm±H cos ( — 



-PmH[i-§:^] (29) 



and can be interpreted as the energy of a magnetic dipole 
with moment Pm = fJ-oT^MsSX. The Zeeman energy of 
the perpendicular magnetization has previously not been 
included in the magnetic energy. It gives a correction to 
the magnetic force on the domain wall. 



dE 
dX 



A A 



hiEx~^x 



Q 



M 



XH 1 



ttA 
2^ 



(30) 



Thus, we include the action of a field component perpen- 
dicular to the wire by replacing the field in Eq. IjlTfl by 
an eflfective field H^q = H {1 + Pm/Qmt). 

We now take into acount the curvature of the wire. 
With decreasing ring radius the angle between neighbor- 
ing spins in the domain wall shrinks. This leads to an ad- 
ditional contribution to the exchange energy of the wall 
when its magnetization points out of the wire plane. 

To calculate the new exchange energy we change the 
spin basis to Cartesian coordinates. To distinguish the 
spin basis from the basis in space we introduce the co- 
ordinates X = COS0, 77 = sin cos ^, and C = sinflsini/). 
Moving along the wire the magnetization performs a ro- 
tation in —9 direction due to the domain wall as well as 
a rotation around the C, axis due to the curvature. For 
small rotations and A/? the Cartesian coordinates are 



X = cos(A/?) cos(6l + M) - sin(A/?) sin(6i + M) cos{4>) 
T] = cos(A/3) sin(6' 4- AO) cos(0) + sin(A/?) cos(6' + AO) 
C sin(6' A6') sin(0). 

(31) 



The exchange energy density is given by 



dx 
dx 



drj 
dx 



dx 



(32) 



From Eq. (EU we obtain 



961 
dx 



-A 
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■ sm 



2 A— — cos ( 

ax r 

2 ^ A 
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(33) 



The first term is equal to the exchange energy density of 
the straight wire. The last term is constant and does not 
depend on the magnetization. In the approximation for 
small (b the other two terms can be rewritten 



AWe 



A—- (2- 

dx r 



^A 



Integration leads to the contribution 



dV AVTcx 



1 AStt 



AS2\\ 



sm 



2AStt 



(34) 



(35) 



of the curvature to the anisotropy energy. The last term 
is a constant which does not depend on X or <j). The 
perpendicular anisotropy energy can be written as 



dVWs.A 



dVK± sin^ e sin^ = K^S2\4? . (36) 



Comparing Eqs. (|35|) and (|36|) one can see that the ad- 
ditional exchange energy due to the curvature can be in- 
cluded into the perpendicular anisotropy by defining an 
effective anisotropy constant 



^icff = ^'^^ + 



A-k 
2A7 



A_ 



(37) 



V. NUMERICAL CALCULATIONS 

To check the applicability of the approximations made 
in our analytical model, i.e. the form invariance of the 
domain wall at small displacements, we have performed 
micromagnetic simulations. We have modelled current 
induced domain-wall oscillations in curved nanowires as 
described in Sec. The current contacts are arranged 
under an angle of 90° to have sufficient distance to the 
domain wall as well as to the ends of the wire (see Fig. 

m 
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We extended the implementation of the Landau- 
Lifshitz-Gilbert-equation in the Object Oriented Micro 
Magnetic Framework (OOMMF)ii by the additional 
current-dependent terms of Eq. Q and implemented 
Runge-Kutta and Adams-Bashforth-Moulton algorithms 
of higher order to speed up the calculations. The cal- 
culations presented here have been performed using the 
explicit embedded Runge-Kutta 5(4) algorithm by Cash 
and Karpi^ The current density is calculated by locally 
solving Ohm's law, thus taking the curvature of the wire 
and the contacts into account. For the spatial discretiza- 
tion a cell size of 1 nm in x- and y-direction and 10 nm 
in z-direction was chosen. Numerical calculations were 
performed for radii of 45 nm, 55 nm, 65 nm, 70 nm, 
85 nm, and 95 nm with different polarized current densi- 
ties jp ~ jP. Wc use the material parameters of permal- 
loy, i.e. the exchange constant A = 13 ■ 10""'^^ J/m and 
the saturation magnetization Mg = 8-10^ A/m. All 
wires have a quadratical cross section S = wt = 100 nm^. 
The applied field in y-direction was chosen to be 125 mT 
to increase the resonance frequency, see Eq. ()20|l. and 
thus to reduce the simulation time necessary for the do- 
main wall to perform several oscillations. Due to the 
small width of the wire this high field has virtually no 
effect on the ground state {H = 0) of the magnetization. 
In the ground state wc obtain a domain-wall width of 
A = 9.25 nm. The difference in the magnetization ori- 
entation 6 between the analytical description of the Neel 
wall and the micromagnetic ground state in the curved 
wire is less than 5°. 

We have determined the eigenmodcs of the magneti- 
zation in the wire by applying a magnetic (5-pulsc in z- 
direction (see Fig. thus exciting all frequencies with 
equal amplitude. To mimic an applied current, the mag- 
netic field pulse has been chosen to point in z-direction so 
that the torque of the field points in the same direction 
as the torque of the applied current [see Eq. (j2J)] . After 
this excitation the system performs damped free oscilla- 
tions. The eigenmodes of the wire are found by spatially 
resolved discrete Fourier transformation (see Fig. (SJi^^i^ 
The higher harmonics and the standing spin waves in the 
wire are neglected in the analytical description. The res- 
onance of the ground mode is observed at a frequency of 

= 15.7 GHz. The higher modes are also indicated in 
Fig. O However in the following wc focus on the ground 
mode. 

We simulated an alternating current with frequencies 
close to the resonance frequency of the domain wall for 
different radii r and Gilbert damping parameters a. Fig- 
ure^ shows the numerically obtained amplitudes for dif- 
ferent radii at fixed a ~ 0.05 and ^ = 0.01. For each 
radius the position and the width of the resonance curve 
have been fitted to the analytical model, Eq. p7|l . to de- 
termine the parameters F{r), uJr{r), and r(r). Note that 
all resonance curves are in excellent agreement with the 
harmonic-oscillator model. The frequencies ujr{r) and 
the damping constants r(r) have been summarized in 
Fig. [SI where they are compared to the analytical cxprcs- 
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00 (GHz) 

FIG. 3: (Color online) Fourier transform M{u}) of the sim- 
ulated magnetization Mx{t) in a curved nanowire with ra- 
dius r = 45 nm and Gilbert damping parameter a — 0.05. 
The wire is excited with a magnetic 5-pulse. The lines show 
the spatially resolved (solid line) and the integral response 
(dashed line). The insets show the spatially resolved discrete 
Fourier transforms for seven selected eigenfrequencies. 
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FIG. 4: (Color online) Amplitude of the domain-wall displace- 
ment versus frequency of the applied current for different radii 
r. The Gilbert damping a = 0.05, the ratio of the exchange 
and spin-flip relaxation time ^ = 0.01, and the polarized cur- 
rent density jp = lO^'^A/m^ are fixed. The crosses denote 
numerical values while the lines are fits with the analytical 
result of Eq. llTTIl . 



sions in Eq. H18|l and Eq. I|20|) . The results coincide if we 
assume Kj_^q = ^^^ + 4^' 7^ ^it^ A'_l = 60000 J/m^ 
for the perpendicular anisotropy [see Eq. H37() ]. The de- 
pendence of the resonance frequency lj^ on the radius r 
according to the phenomcnological model of Saitoh et aliS 
is also shown. It is visible from Fig. [Slthat the analytical 
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1/r (1/nm) 

FIG. 5: (Color online) Resonance frequency ujr and damp- 
ing constant F versus reciprocal ring radius. Shown are the 
values determined from the fits in Fig. |1| (data points) and 
the analytical values (solid lines). The dashed line indicates 
the behavior of the resonance frequency as expected from the 
phenomenological model of Saitoh et ali- 
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FIG. 6: (Color online) Amplitude of the domain-wall dis- 
placement versus frequency of the applied current for different 
Gilbert damping parameters a. The ring radius r — 45 nm, 
the ratio of the exchange and spin-flip relaxation time ^ — 
0.01, and the polarized current density jp = 10^" A/m^ are 
fixed. The crosses denote numerical values while the lines are 
fits with the analytical result of Eq. I17II . 



model and the phenomenological oscillator model yield 
the same eigenfrequencies in the limit of a straight wire 
(r >> 1 /im). For smaller radii the phenomenological 
model gives eigenfrequencies which are significantly lower 
than the ones of the nmiierical calculations. Our analyt- 
ical model including the geometrical corrections fits the 
numerical data very well. 

Figures and [7| show the corresponding data for a 
ring with a radius of 45 nm and different values of the 
Gilbert damping parameter a. The analytical solutions 
are calculated with no free fit parameter. While the data 
points for the damping constant T(a) coincide with the 




FIG. 7: (Color online) Resonance frequency ujr and damping 
constant F versus Gilbert damping parameter a. Shown are 
the values determined from the fits in Fig. |S| (data points) 
and the analytical values (solid lines). 
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FIG. 8: (Color online) Force per wall mass at the resonance 
frequency versus reciprocal ring radius. Shown are the nu- 
merical values (crosses) and the analytical values (line). The 
polarized current density is jp — 10^^ A/m^. 



analytical result, small deviations occur in the resonance 
frequency u>r{a). These deviations can be attributed to 
the finite cell size in our simulations. 

In Figs. |51and|^the values for the fit parameter F{a, r) 
are compared with the analytical result. The analytical 
values exceed the numerically obtained parameters by up 
to a factor of two. This difference has several reasons. In 
Sec. mil we assumed that the ground mode can be de- 
scribed by the motion of the center of the wall X and 
the magnetization angle (p. This neglects spin-wave ex- 
citations and higher wall modes. Calculating the mode 
spectrum excited with a single driving frequency revealed 
a strong coupling between the ground mode and higher 
modes. This coupling is enhanced for small radii. There- 
fore, the force is distributed over several modes, thus de- 
creasing the amplitude of the ground mode. Moreover, in 
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FIG. 9: (Color online) Force per wall mass at the resonance 
frequency versus Gilbert damping parameter. Shown are the 
numerical (crosses) and the analytical values (line). The po- 
larized current density is jp = 10^^ A/m^. 
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FIG. 10; (Color online) Resonance frequency ujr versus 
Gilbert damping parameter a. The data points are the nu- 
merical values obtained for a wire with radius 45 nm and two 
different densities of the polarized current jp. The line is a fit 
according to the analytical result ujr = C/VT+1? from Eq. 
1201 1 with the fit parameter C. 



wires with small radii the current distribution is very in- 
homogeneous with a higher current density at their inner 
edge. We expect that this leads to an additional defor- 
mation of the Neel wall. Another aspect is the finite cell 
size. In the numerical calculations the curved surface has 
been approximated with rectangular prisms. The result- 
ing kinks in the wire wall have a measurable effect on the 
domain-wall motion similar to surface roughness. 



VI. RELATION TO EXPERIMENT 

In the analytical calculations we assume the linear ap- 
proximations sm(X/r) « X/r and sine/) w (j). Nonlin- 



earities cause the deviation of the resonance frequency 
in Fig. IIUI from the analytical form at the current den- 
sity jp = 10^^ A/m^ and Gilbert damping parameters 
below 0.05. Note that these non-linearities are small (< 
0.5%). The current densities at which non-linearities be- 
come important strongly depend on the geometry of the 
wire. Our analytical model allows to derive them for typ- 
ical experimental parameters. In the rest of this section 
we will assume that the damping constant T is not field- 
dependent (^ << J^^j ) and that the squared Gilbert 
damping parameter is small (a^ << 1). These assump- 
tions usually hold in experiments because the domain- 
wall width A is small compared to the radius r and the 
usual values of the damping parameter a are lower than 
0.1. We can express all terms in Eq. (|22|l with the expres- 
sions for r and uJr in Eqs. (|18|l and H20|) when we assume 
that the ratio of exchange and spin-flip relaxation time ^ 
is comparable or less than the Gilbert damping parame- 
ter a. In the case of a non-critically damped oscillation 
(iUr > r) the oscillation becomes nonlinear if the current 
density is approximately 



J = mm 



Tr r^A 



46j 2aujrbj 



(38) 



The experimental current densities which Saitoh et aliS 
applied on a wire with cross section S = 3150 nm^ and 
radius r = 50 /xm are well below this current density. 
They have determined a domain-wall width A = 70 nm, 
a domain-wall mass m = (6.55 ± 0.06) • 10"^"^ kg, and a 
domain-wall relaxation time t = ^ = (1.4±0.2) • 10~^ s. 
To compare these findings with our analytical model we 
have to make an assumption for the effective anisotropy 
constant Kj_qQ- Assuming the value K^^^q ~ 60000^^-1- 
om the fit in Sec. Eq. l(T5|) leads to the 
1.55 • 10"^'^ kg. This mass is comparable to 
the experimental value, despite the uncertainty of Kj^^q 
due to the different wire dimensions. 

As mentioned in Sec. IIIII the analytical calculations 
lead to relations between the micromagnetic material 
parameters and the parameters of the harmonic oscil- 
lator. These can be used to experimentally determine 
the Gilbert damping parameter a from the experimental 
data. From Eqs. (|18|l . (|20|l . and H15|) one can derive the 
relation 



A 

2Ar r- 

mass m 



2T-fHX mTjX 



HoMsS 



(39) 



With the domain-wall mass and the domain-wall relax- 
ation time of Saitoh's experiment we get a Gilbert damp- 
ing parameter of a = 0.0114 ± 0.0017. This value agrees 
quite well with the experimental values of Nibarger— 
and Schneider— which range from 0.008 to 0.017 for film 
thicknesses between 10 nm and 93 nm. 
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VII. CONCLUSION 

The current-induced motion of a domain wall in thin 
curved nanowires has been investigated. A harmonic- 
oscillator model which so far had only been introduced 
phenomenologically is derived from the LLG equations 
extended by the spin torque according to Zhang and Lij2i 
This derivation relates micromagnetic material parame- 
ters to the characteristic quantities describing the oscil- 
lating domain wall under the influence of an alternating 
driving current. It is shown that the dipole moment of 
the wall as well as the curvature of the wire have an im- 
portant influence on the resonance frequency and damp- 
ing constant of the oscillation. The domain wall can be 
seen as a quasi particle in a parabolic potential well which 
is acted upon by a current-induced force. The phase and 
magnitude of the force depend on the frequency of the 
current. The analytical results have been compared to 
numerical simulations. They agree very well. Our analyt- 
ical solution suggests new methods to determine material 



parameters which are otherwise difficult to measure, e.g., 
the non-adiabatic term of the spin torque can be deter- 
mined from the phase shift between the applied current 
and the overall magnetization. Moreover, the Gilbert 
damping parameter a and the domain-wall mass m fol- 
low from a measurement of the resonance frequency LUr 
and the damping constant T of the oscillations. 
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